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Abstract 

Exit times for stochastic Ginzburg-Landau classical field theories with 
two or more coupled classical fields depend on the interval length on which 
the fields are defined, the potential in which the fields deterministically 
evolve, and the relative stiffness of the fields themselves. The latter is 
of particular importance in that physical applications will generally re- 
quire different relative stiffnesses, but the effect of varying field stiffnesses 
has not heretofore been studied. In this paper, we explore the complete 
phase diagram of escape times as they depend on the various problem 
parameters. In addition to finding a transition in escape rates as the rel- 
ative stiffness varies, we also observe a critical slowing down of the string 
method algorithm as criticality is approached. 



1 Introduction 

In a previous paper [T] (hereafter GS), the authors introduced and solved a 
system of two coupled nonlinear stochastic partial differential equations. Such 
equations are useful for modelling noise- induced activation processes of spatially 
varying systems with multiple basins of attraction. Examples of such processes 
include micromagnetic domain reversal [HG2, pattern nuclcation 0HS], transi- 
tions in hydrogen-bonded ferroelectrics [TJ, dislocation motion across Peicrls 
barriers [5], and structural transitions in monovalent metallic nanowires [UlllOj. 
It is the last problem in particular that the model introduced in GS was con- 
structed to analyze. 

The GS model provided a mathematical realization of a stochastic Ginzburg- 
Landau field theory consisting of two coupled classical fields, denoted 4>\{z) and 
4>2(z), defined on a linear domain of finite extent L. Stochastic partial differ- 
ential equations of this type are constructed to model noise-driven transitions 
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between locally stable states. In the especially important case of weak noise, 
where the transition rate is of the Arrhenius form Toe^ AE ^, with prefactor 
Tq and activation barrier AE independent of the noise e, the transition path 
occurs near (i.e., within a lengthscale of order 0(e 1 ^ 2 )) the saddle (or col) of 
least action connecting the two stable states. 

The two-field model displayed several interesting features, including a type 
of "phase transition" in activation behavior as L varied. The transition was 
driven by a change in the saddle state, from a uniform configuration at small L 
to a spatially varying one ( "instanton" ) at larger L. This transition has been 
noticed and analyzed for the case of a single field [TTT[T?| , but had not been seen 
in the rarely studied case of a system with two coupled fields. Perhaps more 
remarkably, the system admitted an exact solution for the instanton state; such 
exact solutions are rare in the case of nonlinear field theories with a single field, 
much less a nontrivial system of coupled fields. 

The introduction of two fields was required to study transitions among dif- 
ferent quantized conductance states in non-axisymmetric nanowires. The ax- 
isymmetric case had previously been treated theoretically in [5J[T3]. However, 
detailed studies using linear stability analysis by Urban et al. |14j indicated that 
roughly 1/4 of all such transitions involved cither non-axisymmetric initial or 
final states, or else a least-action transition passing through a non-axisymmctric 
saddle. To describe such transitions, one field (<j>i(z)) describes radial variations 
along the wire length and the other {<j) 2 {z)) describes deviations from axisym- 
mctry. 

One restriction of the analysis in GS was that the respective bending coef- 
ficients «i and ft2 of the two fields were taken to be equal. However, this is 
generally not the case in real nanowires [14] . Therefore, in order to apply the 
model to actual transitions, as well as to provide a complete picture of the ac- 
tivation behavior in such systems, we need to consider the case where K\ =/= k?. 
In such cases analytical solutions cannot be found and we need to rely on nu- 
merical methods. The study of this more general problem is the subject of this 
paper. 
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2 The Model 



Consider two coupled classical fields 4>i(z), 4>2(z) on the interval [— L/2,L/2], 
subject to the energy functional 

L /2 

h= f (^iM(f+k^w) 2 +%,w)^. (i) 

J-L/2 Z Z 

where 

J7(0i,0 2 ) = -y^i + J&i - y02 + ^02 + ^102 ( 2 ) 

The bending coefficients K\, k 2 can be related to the wire surface ten- 
sion [51113] . The arbitrary positive constants fj,±, [i2 arc chosen such that [i\ ^ fj,2, 
breaking rotational symmetry. (The case [i\ = /i 2 has been investigated ana- 
lytically by Tarlie et al. [15] in the context of phase slippage in conventional 
superconductors. ) 

If the system is subject to additive spatiotemporal white noise, then their 
time evolution is governed by the pair of stochastic partial differential equations: 

4>i = K 1( f>'{ + /!i</>i - - 4> 1( f>l + V2e^(z, t) 

4>2 = K2<t>2 + ^2^2 - 4>2~ <t>\ 02 + V2e&(z,t) , (3) 

where £i j2 are the spatiotemporal noise terms satisfying < t\)^j{z2, i 2 ) >= 

8{z\ — Z2)8{t\ — t2)5ij 1 i,j = 1,2. If the noise is due to thermal fluctuations, 
then by the fluctuation-dissipation theorem e = ksT . 

The activation energy AE and prefactor T in the Arrhenius rate formula 
depend not only on the details of the potential ([2]), but also on the interval 
length L on which the fields are defined, and on the choice of boundary con- 
ditions at the endpoints z = — L/2 and z = L/2. It was shown in |16) that 
Neumann boundary conditions are appropriate for the nanowire problem, and 
we will employ them here. 

In the usual case of a single field, the bending coefficient k plays a role 
in determining the intrinsic lengthscale (and therefore the transition length at 
which the saddle state changes) on which field variations occur; but once it is 
absorbed into a dimcnsionless lengthscale by rescaling along with the parameters 
determining the potentials, it plays no further role. Now, however, there are 
two bending coefficients, and varying their relative magnitudes can in principle 
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lead to new phenomena. The aim of this paper is to study the effects of these 
variations. 

The metastable and saddle states are time-independent solutions of the zero- 
noise equations: 

Kx<t>'l = "Mi 01 +<Pi+ <fii4>l 

= -M202 + 0| + fa<t>\ (4) 

Without loss of generality, we choose fa > fa. 

Then there are two metastable states: (f>i fS = iy/ToT, 4>2, s = 0; two spatially 
uniform saddle states: </>i jU = 0, 02, u = ^y/fa] arL d spatially nonuniform saddle 
states, or instantons. When k\ = K2{= 1), analytical solutions for the instanton 
saddle states can be found: 



4>\%,{z) = ±V™V (2^1 - A*2) - tMjjli - fa)sn(yjfa - fa z\m) (5) 
<A2 n m( z ) = ±VM2 - m (Mi - M2)dn( V/ii - M2 (6) 

where sn(.|m)and dn(.|m) are the Jacobi elliptic functions with parameter m, 
whose periods are AK(m) and 2K(m) respectively, with K(m) the complete 
elliptic integral of the first kind [T7]. The parameter m G [0, 1] is related to 
interval length L through the Neumann boundary conditions, with m — > + 
corresponding to L — > L+, where L c is the critical length, and m — > 1 corre- 
sponding to L — > oo EHIHJQj]- When K\ = K2 = 1, 

L= 2K ^ (7) 

Wc found in GS that varying L triggers a transition between the uniform 
and instanton saddle states; the critical length L c is determined by: 

VMi - M2 VMi - M2 

This results in a transition in the activation behavior, including anomalous 
behavior at the critical length. Such a transition may have already been seen 
experimentally |10j , in a crossover from ohmic to nonohmic conductance as the 
voltage across gold quantum point contacts increases |18| . We will show below 
that the same effect occurs when the ratio Ki/«2 is varied. 

As noted above, the transition rate in the low-noise (e — > 0) limit is given 
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by the Kramers formula: 



r-T cxp(-A£/e) . 



(9) 



Here AE is the activation barrier, that is, the difference in energy between the 
saddle and the starting metastable states, while To is the rate prefactor: 



In the above equation A s is the linearized dynamical operator describing pertur- 
bations about the metastable state; similarly A u describes perturbations about 
the saddle. A u .i is the single negative eigenvalue of A u , corresponding to the 
direction along which the most probable transition path approaches the saddle 
state. The behavior of Tq becomes anomalous at the critical point L c , where 
fluctuations around the most probable path become large. 

3 Calculation of the Minimum Energy Path 

Computation of exit behavior requires knowledge of the transition path(s), in 
particular behavior near the local minimum and the saddle. In our model, both 
are found as solutions of two coupled nonlinear differential equations pQ. A 
powerful numerical technique constructed explicitly for this type of problem is 
the "string method" of E, Ren, and Vanden-Eijnden [THEO]. The algorithm 
proceeds by evolving smooth curves, or strings, under the zero-noise dynamics. 
These strings connect the beginning and final locally stable states, and in be- 
tween the two ends each string contains a series of intermediate states called 
"images". The method is constructed so that the string evolves towards the 
most probable transition path. The evolution proceeds until the condition for 
equilibrium is reached: 



where H is given by ([!} and [5'H] 1 ~ is its component perpendicular to the string. 

Once equilibrium is reached, the string images correspond to the configura- 
tions sampled by the system at different stages of the activation process. The 
image with highest energy is the one nearest the saddle state. In order to get an 
accurate result, the distribution of images needs to be sufficiently fine-grained. 
In our computation, we used 61 images (including the two ends of the string); 




(10) 



[snr = o 



(ii) 
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because of the symmetry of our energy functional, the image in the middle 
corresponds to the saddle. 

When such symmetry is absent and the location of the saddle needs to be 
determined with high precision, one can use an alternative method to the brute 
force one of simply increasing the number of images along the string. This 
alternative requires first finding a rough approximation of the most probable 
transition path, again using the string method but with a small number of 
images, and then switching to a "climbing image" algorithm in which one picks 
up an image that is believed close to the saddle and then drive it towards the 
saddle. The climbing force is obtained from inverting the energy gradient along 
the direction of the unstable eigenvector of the saddle state. Details can be 
found in Q33I2S]- 

Wc have found an analogue to critical slowing down in the current context: 
near criticality convergence of the string method becomes increasingly slow. 
Expanding the energy functional around the saddle reveals that the lowest stable 
eigenvalue vanishes to second order, leading to a rapid increase in relaxation 
time. This phenomenon will be further investigated in the following sections. 

4 Results 

We now turn to the case ki ^ K2. To begin, we fix K2 = 1 and vary k\. We 
consider the cases where K\ is both less than and greater than 1. Because the 
critical length now depends on K\, we denote it L c (ki). 

As noted earlier (cf. ([8]), L c (l) = -^=== : below L c (\) the saddle is spatially 
uniform, and above L c (l) it is spatially varying pQ. The situation becomes more 
complicated when K\ ^ 1. Fig. [T] summarizes our results when fi± = 3, fi2 = 2 
and L > L c (l) = it. In this and Fig. [3J the saddle state (whether uniform or 
instanton) is denoted </>f addlc (z), i — 1, 2. Wc find that as k% increases, the spatial 
variation of the instanton becomes increasingly suppressed until the instanton 
finally collapses to the uniform state. Conversely, when L < L c (l), the instanton 
state is retrieved for n\ < 1 (cf. Fig. [2]). 

This effect can be understood as follows. In the limit of low noise, the 
transition occurs over the saddle state of least energy. An increase in k\ raises 
the bending energy of any nonuniform configuration, and therefore that of the 
instanton, while leaving the energy of the uniform saddle unchanged. When 
the energies of these two states cross, the activation process switches from one 
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saddle state to the other. This is seen explicitly in Fig. [31 where we plot the 
energy of the saddle state against K\ for both L = 0.25 and L = 4.51. In 
these figures the curve to the left of the dashed line is the instanton branch, 
which increases monotonically until it reaches a constant value: the energy of 
the uniform saddle state. 

We next investigated the question of whether the transition from uniform 
saddle to instanton (or vice- versa) as Ki varies occurs as a continuous crossover 
or as an abrupt phase transition. If the latter, then we also need to determine 
the order of the transition. 

In pp, the uniform — > instanton saddle transition was induced by changing 
L at fixed n\. There we concluded that the transition was reminiscent of a 
second-order phase transition, in the asymptotic e — > limit. This follows from 
the continuity of the activation energy at all values of L, including L c (l). (For 
examples of potentials where the activation energy jumps at a precise value of 
L, corresponding to a first-order transition, see [21].) In fact, it can be shown 
analytically that the first derivative of the activation energy curve with respect 
to L is also continuous everywhere, but the second derivative is discontinuous 
at L c {\). 

Similarly, Fig. [3] suggests that there is indeed a continuous phase transition, 
in that the activation energy changes continuously as one passes through the 
transition, as ki varies for fixed L. This continuity leads to a divergence in the 
transition rate prefactor, shown in Fig. [4] (similar to that induced by changing 
L at fixed K\ in Q]). The value of ki c (L = 4.51) where the prefactor diverges 
and that where the energies of the respective saddles cross agree to within a 
numerical error of 10~ 2 . 

What causes this divergence? Away from criticality, the spectrum of the 
linearized dynamical operator A„ about the saddle consists of a single negative 
eigenvalue, whose corresponding eigenvector determines the unstable direction, 
with all other eigenvalues positive. As criticality is approached, the smallest 
positive eigenvalue, denoted A Uj 2, approaches zero. This signals the mathemat- 
ical divergence on the "normal" lengthscale of 0(e 1 ^ 2 ) of fluctuations about 
the saddle, and by ([T0| is seen to lead to divergence of the prefactor. (For a 
discussion of how to interpret this "divergence", see |2"2"j.) 

The eigenvalue spectrum about the uniform saddle can be analytically cal- 
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culated [T]. The eigenvalue A Uj 2 is found to be 



A u ,2 = -jj- - (in - Hz) ■ (12) 

At fixed L, this switches from negative (unstable) to positive (stable) as Ki 
increases, as shown in Fig. [5] This change of sign corresponds to a transition 
from an instanton saddle to a uniform one as K\ varies. Using this approach, 
the curve L c vs. ki c can be derived analytically as the locus of points where 
A„,2 — and thus the full phase diagram determined as represented by the solid 
curve in Fig. [6l 

Wc have also studied the behavior of the transition rate prefactor in a wide 
range of values of L numerically, all of which lead to the same conclusion as 
described above. Fig. [7] shows the divergence of To at different L c 's and their 
corresponding Ki c 's. 

5 Discussion 

We have solved the general two-field model of (U} and © for its full parameter 
space. We have uncovered a new mechanism for the transition in the switching 
rate, and shown that it has features of a second-order phase transition. 

In the one-field case, the mechanism behind the transition is not difficult to 
understand. At smaller L (recall that this is a dimensionless lengthscale, in units 
of a reduced length that includes the single bending coefficient k), bending costs 
(in conformity with the boundary condition) are prohibitive, and the uniform 
saddle therefore has lower energy than the instanton. At large lcngthscales, 
the uniform saddle has a prohibitive bulk energy (i.e., potential difference with 
the stable state), whereas the instanton differs from one or the other stable 
state only within the domain wall region, whose lengthscale remains of 0(1). 
What is perhaps less intuitive is that the transition in saddle states should be 
asymptotically (as e — > 0) sharp. 

Here wc have uncovered a second mechanism for the transition to occur: as 
noted in Sec. SJ increase of K\ when L > L c (l) suppresses spatial variation, 
causing the instanton (again in a sharp transition) to collapse to the uniform 
saddle. Conversely, the instanton state can be retrieved for L < L c (l) when n\ 
decreases] of course, bending becomes increasingly favorable energetically. The 
result is a phase diagram in (L, ki/k 2 ) space, as in Fig. [SJ where a phase sepa- 
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ration line denotes the boundary between the uniform and instanton "phases" . 

We close with some remarks about the string method as applied to this 
problem. 

A randomly placed string will relax towards the most probable transition 
path along the stable direction of the saddle. In Sec. @] we defined the small- 
est positive eigenvalue (corresponding to the stable direction) of the linearized 
operator A u . As a second-order phase transition is approached, X u ^ drops to 
+ , so that the energy landscape curvature in the stable direction becomes very 
small. When the string arrives in its neighborhood, the restoring force exerted 
along its normal direction correspondingly becomes small leading to slow con- 
vergence. If one sits right at the critical point, the string will not arrive at the 
saddle. 

The string method assumes that most of the probability flux from the re- 
actant to the product state is carried by one (or more generally a few) paths 
through the saddle state, in each of which the probability flux is tightly confined 
to a narrow quasi - ID region in state space. However, near criticality the path 
splays out in the direction normal to the longitudinal transition path. In this 
case one needs to consider transition "tubes" , inside which most of probability 
flux is concentrated. The equilibrium condition (|3|) corresponds to conditions 
away from criticality, where the transition tube is thin. 

The equation for the path of maximum flux is derived in |23| . where it 
is noted that the reaction flux intensity must be maximized along the thin 
transition tube (or the string, when using the string method). An alternative 
derivation can be found in |24j . 
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K\ = 2.20 



= 2.20 
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(a) 0f ddle (z) 



(b) 0| addle (z) 



Figure 1: The saddle states (j)*^^ passed by the system at different relative field 
stiffnesses (= K\). The evolution of the saddle is described as the continuous 
transition of colors. K\ was increased from 1.95 to 2.20 with an increment of 0.01. 
The arrows indicate the suppression of the instanton at L > L c (l) = f as 
K\ increases. Here fix = 3, \i-2 = 2, and L = 4.51. 





(a) </>f ddle (z) 



(b) 0| addle (2) 



Figure 2: The saddle states ^fi^Jf passed by the system at different relative 
field stiffnesses(=Ki). The evolution of the saddle is described as the continuous 
transition of colors. k,\ was decreased from 1.0 to 0.5 with a decrement of -0.01. 
The arrows indicate the retrieval of the instanton at L < LAV] = , £ „ as Ki 

Ly ' v/>l-/i2 

decreases. Here H\ = 3, ^2 = 2, and L = 0.25. 
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0.2 0.4 0.6 0.8 1.0 1.95 2.00 2.05 2.10 2.15 2.20 



(a) L = 2.50 (b) L = 4.51 

Figure 3: Energy of the saddle state as a function of «i, where = E sa ddie — 

^metastable ■ 

The dashed lines indicate that n lc s=* 0.63 for L = 2.50 and k 1c sa 
2.06 for L = 4.51, where energies of the instanton and uniform saddles cross. 
The region k\ < k\ c (L) corresponds to the instanton saddle, and K\ > Ki c (L) 
to the uniform saddle. Here [i\ = 3, and [l% = 2. 
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Figure 4: Behavior of the prefactor Tq as calculated numerically using (ITU1) when 
Hi = 3, /12 = 2, and L = 4.51. 
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Figure 5: Smallest nonnegative eigenvalue A Uj 2 (solid line) of the saddle state 
for iii = 3, fj,2 = 2 and L = 4.51. For ki < 2.06 this eigenvalue corresponds to 
the instanton saddle, and for K\ > 2.06 to the uniform saddle. The extended 
dashed line shows the continuation of this eigenvalue for the uniform saddle in 
its unstable regime below k\ c . 
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2.5 3.0 3.5 4.0 4.5 5.0 



Figure 6: The phase diagram at /ii = 3, \x<z = 2. The dots represent numerically 
determined values where the energies of the uniform and instanton saddles cross. 
The solid line was computed analytically using (fT2"|) . that is, by determining the 
relation between k\ and L along which the smallest nonnegative eigenvalue of 
the uniform saddle is zero. 
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Figure 7: The divergence of Tq at L ranging from 2.6 to 5.0. Here n\ = 3, 
(j,2 = 2, and the values of K\ c (or the corresponding L c ) correspond to the dots 
in Fig. El The peaks are of different heights because the speeds of divergence 
are not necessarily the same for every L. 
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